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Abstract 

Given asymptoticahy flat initial data on M'^ for the vacuum Einstein field equation, 
and given a bounded domain in Af, we construct solutions of the vacuum constraint 
equations which agree with the original data inside the given domain, and are identical 
to that of a suitable Kerr slice (or identical to a member of some other admissible family 
of solutions) outside a large ball in a given end. The data for which this construction 
works is shown to be dense in an appropriate topology on the space of asymptotically 
flat solutions of the vacuum constraints. This construction generalizes work in 
where the time-symmetric case was studied. 

1 Introduction. 

We study the asymptotics of solutions to the vacuum constraint equations. In particular, 
we show how to approximate given initial data for the Einstein vacuum equation by data 
which is exactly that of a spacehke sUce of a suitably chosen Kerr metric outside a compact 
set. For the time-symmetric case, the constraints reduce to the equation R{g) = 0, and it 
was shown in [HI how a Schwarzschild exterior can be selected and glued to asymptotically 
fiat (AF) time-symmetric data and compactly perturbed to a solution of the time-symmetric 
constraint. In the non-time-symmetric case, one has to consider the second fundamental 
form, and hence linear and angular momentum of the data at infinity, for which the Kerr 
solution (for example) is aptly suited. We note that these constructions show that unique 
continuation fails in a serious way, but it also points to the reahty that a compact piece of 
the data determines (in some sense) only part of the asymptotic structure, and indeed that 
which it determines is only the part that is governed by the energy-momentum and angular 
momentum. Having vacuum initial data which is identical to Kerr data outside a compact 
set is important for understanding the global evolution for the Einstein equations. In partic- 
ular, such data evolves to produce a spacetime with particularly nice behavior at null infinity. 

We briefiy recall the basic setup of the constraint equations. A solution of the vacuum 
Einstein equation is a Lorentzian manifold {S'^,g) satisfying 

mcig)-^R{g)g = 0. 

In this equation, Ric{g) and R{g) denote the Ricci and scalar curvatures, respectively, of the 
metric g. By taking the trace we see that this is equivalent to Ric{g) = 0. It is well-known 
that Einstein's equation admits an initial value formulation, in which the vacuum initial 
data consist of an oriented three-manifold M^, a Riemannian metric g and a symmetric 
(0, 2)-tensor K on M. The Gauss and Codazzi equations provide the constraints upon g and 



K in order that they form, respectively, the induced metric and second fundamental form of 
M inside a Ricci-fiat spacetime {S, g) ; in the vacuum case these constraints are |Wj 

R{g) - IKI"^ + H"^ = 
diVg{K)-VH = 0. 

Here H = Trg{K) = g^^Kij denotes the mean curvature, and all quantities are computed 
with respect to g. 

We rewrite these equations by introducing the momentum tensor 

^ij = K'^ -Trg{K)g'K 
We also introduce functions Ti and $ by 

H{g,Ti) = R{g) + ^{Trgnf-\n\' 
*(5',7r) = {n{g,iT),diVgTT). 

The constraints then take the form $((7, tt) = 0. 



Two well-known solutions of the vacuum Einstein equation are the Schwarzschild solution 
and the Kerr solution. The Schwarzschild metric is characterized by rotational symmetry 
and is static; indeed outside the horizon we have coordinates in which the Schwarzschild 
metric takes the form 

gs = -(l-'-^^ dt- + (1 - V) " + 

In particular we see that the metric is static, timelike Killing field which is orthog- 

onal to the time-slices. The Kerr solution is axisymmetric and stationary; the analogous 
Killing field is not orthogonal anymore to the time-slices (physically, the black hole is rotat- 
ing), and this can be easily seen in Boyer-Lindquist coordinates |MTWj . |Wj . 

What is more important for us is that these solutions are actually families of solutions 
of the Einstein equations. In fact by thinking of a fixed asymptotically fiat coordinate 
system near infinity in the spacetime, by varying the total mass and angular momentum, 
and considering families of spacelike slices of these metrics, we get a ten-parameter family 
of asymptotically fiat solutions of the vacuum constraint equations near infinity in M?. The 
purpose of considering this family is to control the energy-momentum [E, P) and the angular 
momentum J at (spacelike) infinity, as well as a quantity C related to the center-of-mass. In 
an asymptotically flat chart these quantities are deflned as limits of integrals over euclidean 
spheres with surface measure dag and outer normal v taken with respect to g 

^ = "i4z y i9ij,i - 9ii,j) T^^dag (1) 
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Pj = — lim (v "S^ Tx^^v^ da„ (2) 
C*^ = lim i x'' {g^j,i- g^iJ)u^ dag - Mm i y^ig^ky'-gii^^^)dag (4) 

_R— >oo J ' ' _R— >oo J ' ' 

|x|=i? ^'-^ |x|=_R * 

In the case when g is AF, conformally flat near infinity and has zero scalar curvature, the 
quantity C is proportional to mco, where m is the mass and Cq is the coordinate translation 
which makes the |x|~^-terms in the expansion of the conformal factor vanish [Cj. We note 
that it may be useful to shift the quantity C (for example an initial shift of coordinates will 
do this), and also we note that the vector fields Yi are the basic rotation fields in R^, for 
example 

Yi = a; 7— - x^ 



dx^ dx"^ 

The integrand for the angular momentum is essentially the cross product of the position 
vector with the momentum tensor, in analogy with the angular momentum of classical me- 
chanics. 



We can prescribe conditions under which the above asymptotic integrals converge for 
solutions of the vacuum constraints. These conditions are asymptotic even/odd conditions 
on the metric tensor and second fundamental form of the initial data; we call this condition 
(AC). The conditions are 

(7,,(x) = 5,, + 0(1x1-1) ^ 0(|x|-2) 

<7.,(x) - g,,{-^) = Odxr^) ir,,(x) + i^,,(-x) = 0(1x1-3) 

5(ij-fc(x) + fc(-x) = 0(|x|-3) Kij^ki^) - Kij-fe(-x) = 0(|x|-^) 

We require analogous conditions on successive derivatives as needed. It is also worthwhile 
to remark that under these asymptotic conditions on solutions of the constraints, the corre- 
sponding boundary integrals above can also be computed with respect the euclidean metric, 
and the limiting values are the same, a fact we may use without further comment. 

In Section 3 of this paper we show that the vacuum initial data sets which satisfy (AC) 
are dense in a suitable weighted Sobolev space {gij — 6ij G W'^f and vr G W^f_g for p > 3, 
6 G (1/2, 1)) in the set of all vacuum initial data which satisfy the standard decay assump- 
tions (but not the asymptotic symmetry). This space is strong enough that the total energy 
and linear momentum are continuous on the space. (Of course, the angular momentum is not 
continuous, nor does it appear to be well-defined on this more general space.) This theorem 
is an analogue for the full constraint equations of the theorem of Schoen-Yau |SYj in the 
time-symmetric case {K = 0). It is however much more subtle to prove because the general 
constraint operator is not as well understood as the time-symmetric operator. We show that 
it is possible to achieve a very special type of asymptotic behavior which in particular satis- 
fies (AC) . These special asymptotic conditions require that outside a compact set there exist 
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a positive function u and a vector field X so that Qij = u'^Sij and = u'^{LxS — divs{X)Sij) 
where LxS denotes the Lie derivative of the euchdean metric 6 with respect to X. It is not 
difficult to see that such data automatically satisfies (AC). An important feature of these 
asymptotic conditions is that the total energy and momenta (both linear and angular) may 
be read off the asymptotics, and the behavior of these quantities directly effects the asymp- 
totic geometry. This type of asymptotic behavior will be discussed in detail by the second 
author in [5]. The proof of the density of solutions involves taking an arbitrary solution of 
the vacuum constraint equations and patching g to the euclidean metric and tt to in a 
large annulus. It is then necessary to reimpose the constraints by solving a PDE for {u,X). 
The cut-off data may be viewed as an approximate solution, and one may attempt to use 
the inverse function theorem to correct it to a solution. The problem which arises is that 
the corresponding linearization is not clearly an isomorphism in the natural spaces. This 
problem is overcome by allowing the extra flexibility of addition of a suitably chosen family 
of compactly supported deformations (of (y',7r). In order to make this work, we need to use 
the fact that the constraint operator (at an arbitrary asymptotically flat initial data) has 
surjective linearization in these weighted Sobolev spaces. This result is proven for maximal 
data by Choquet-Bruhat, Fischer, and Marsden in |CFMj . and the general case is due to 
Beig and O Murchadha |BOj . 

We now outline the basic approach to the main theorem. Given asymptotically flat ini- 
tial data on M satisfying (AC) in an appropriate chart at infinity in a given end, we take a 
"large" radius R and within the annulus from R to 2R, we smoothly patch the given metric 
and second fundamental form to the metric and second fundamental form coming from a 
slice in Kerr, or from another suitable family satisfying (AC) near infinity. This will produce 
an approximate solution of the vacuum constraints. The approximate solution is altered 
with a smooth perturbation compactly supported within the closed annular gluing region 
to data {g,Tf), whose constraint function ^{g,7i) lies in a finite-dimensional vector space of 
dimension ten; of course we want ^{g,7i) to be zero. (Note that outside the annulus the 
constraint function vanishes by design.) We then show that there is a choice of data to 
glue on so that the resulting perturbed metric and second fundamental form are solutions 
to the constraint equations. The ten-dimensional obstruction space is countered by the ten 
degrees of freedom afforded by {E, P, J, C); in fact the method proves that suitable models 
near infinity are precisely those for which we can effectively attain values of {E, P, J, C) by 
moving within the family. As in the previous paper, the key to the proof is exploiting the 
overdetermined-ellipticity of the adjoint of the linearization of we include a fairly com- 
plete discussion of the analysis, but some simple details are quoted from 0. The conclusion 
of the proof is the observation that the map from the parameter space to the obstruction 
space, a map between two ten-dimensional Euclidean spaces, has non-zero degree; this map 
is more complicated than before (where only four parameters were needed), but recognizing 
the (AC) condition affords some economy in computing this map as compared to [C], where 
an expansion of the conformal factor is used. 

We remark that the results of this paper are easily extended to allow M to have multiple 
asymptotically flat ends. For convenience of notation we write the proofs under the assump- 
tion that M has only one end. 
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The results of this paper (with the exception of Section 3) were announced in the spring 
of 2000, and have been widely communicated to the mathematical GR community. In the 
meantime, P.T. Chrusciel and E. Delay |(]Dj have obtained a version of these results as 
well. Their paper employs our basic techniques (and those of the first author in [HI) to give 
a number of interesting applications; in addition their paper gives an elegant and explicit 
description of the 10-parameter Kerr family of initial data which we do not do here. 

2 Preliminaries. 
2.1 Basic Notation. 

Let Q denote a compactly contained domain in a smooth three- manifold M^. Unless noted, 
we assume the boundary is smooth. We list here some notation, and we define some function 
spaces which we will find useful. 

• Ric((y') = Rij and R{g) = g^^Rij denote the Ricci and scalar curvatures, respectively, of 
a Riemannian metric g on M. We use the Einstein summation convention throughout, 
as well as the convention of using a semicolon to denote covariant differentiation and 
a comma to denote partial differentiation, and our convention for the Laplacian is 

• We let diig denote the volume measure induced by g, dag the induced surface measure 
on submanifolds, dx the Lebesgue measure on Euclidean space, and d^ the Euclidean 
surface measure. 

• We denote by the Hilbert space of tensor fields which are square-integrable along 
with the first k weak covariant derivatives, with the standard if'^-inner product induced 
by the metric g. It will be clear in context which type of fields are being discussed, if;^^ 
denotes spaces of tensors which are in H'' on each compact subset. We may abbreviate 
products with mixed regularity by using superscripts, like H^'^ = H'^ x H^. 

Similarly we define the Holder spaces C*^'". 

• We may wish to use symbols to distinguish various types of tensors. For example, let 
AA^ > |) denote the open subset of of Riemannian metrics, and J\A^'°' denotes 
the open subset of metrics in C'^'". iS(2,o) denotes the space of symmetric (2,0)-tensor 
fields, and X denotes the space of vector fields. We use superscripts as above to denote 
the desired Sobolev or Holder regularity of the fields. 

• Let p be a smooth positive function on fl. Define L^p{^) to be the set of locally L'^{dfig) 
functions / such that fp^^"^ G L'^{Vl,dpg). The pairing 

makes L'^p{^) a Hilbert space. It is clear how to extend this definition to higher order 
tensors. 
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• Let Hp{Q) be the Hilbert space of tensor fields in L'^{^) along with the first k covari- 
ant derivatives, the inner product defined by incorporating the L^(n)-pairings on the 
covariant derivatives of order 0, . . . ,k. 

• We now want to define the weighted Holder spaces C^l1{Q) (0 < a < 1). We will 
consider p which near dfl decays as a power of or exponentially in the distance to the 
boundary. The weighted Holder space is defined as the subspace of C'^'"(f2) comprised 
of functions / for which the norm 

||/|U,a,p-i •= \\fP~^\k,a 

is finite; this is a Banach space. (Unless noted otherwise, norms will be taken over Q.) 
The following lemma concerning the density of H'^^fl) in Hp{Q) will be useful. 

Lemma 2.1 Assume that p is bounded from above. For k > 1, the subspace H'^{Q) (and 
hence C°°{^1)) is dense in H^{Vi). 

Proof: For small r > let VL^ be a slightly expanded domain compactly containing Vl. Let 
: fi^ — >• f2 be a diffeomorphism which is C'^-close to identity map. Now suppose that / G 
Hp{n) and let = foFr. Note that the restriction of to f2 is in H''(fl). We will show that 
fr^fim Hp{Q). First note that by the chain rule we have |-D"/t| < (^J2\i3\<a l-^^/l for 
I a I < k. Now, the dominated convergence theorem implies that lim Jq^^q^ \D°'f\'^pdpg = 

for |a| < A; where Qs denotes the subset of Q consisting of points at least a distance 6 from 
the boundary. It follows that for r > we also have lim L, „ \D°'fJ'^pda„ = 0. Thus, given 

e > we may choose 5 so small that (/j^^j^^ \D"' fr — D°' f\^ pd^gY^'^ < e/2, and then choose 
T SO small that (J^^ |-D"/t — D"' f\'^ pdpgY^"^ < e/2. Summing these completes the proof. ■ 

2.2 Linearization of the Constraint Map. 

We gather here a few facts we will use in what follows |FMlj . |FM2j . 

Lemma 2.2 The scalar curvature map is a smooth Banach space map, as a map R : 
^ H\Q) (I + 2 > I), or R : M'^+^'^iU) C'^'^ill) {k > 0). The lineariza- 
tion Lg of the scalar curvature operator is given by 

Lg{h) = —Ag{Trg{h)) + divg{divg{h)) — h ■ Ric{g) 

in either of the above spaces. The formal L"^ -adjoint L*g of Lg is given by 

L*g{f) = -{Agf)g + Hessgif) - f Ric{g). (6) 
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Lemma 2.3 The constraint map $ is smooth Banach space map, as a map $ : x 

^Imi^) ^ X (i + 2 > I), or ^ ■. TW^^+^.ap ^ ^JJfy^p ^ c^'-^p x 

(k > 0). The formal L"^ -adjoint D^*^^^-^ of the linearization is given by 

D%^^){f,X) = DHl^^if) + Ddtvl^^iX) where 

D-Hl^Af) = ((^*/)., + ((Tr,7r)7r,,-27r,,7r5)/, {{Tr,n) g^^ - 27r^^) f) 
Ddivl^^^^ (X) = 1 ((Lxvr),, + (Xj)7r,, - (X,7r|^, + vr^.X, ) 
-Xk-mT^^"'gij - Xk-K^^gij, -{LxgY^^ . 

We point out here that the key terms we will use are the L*/-term in the first component 
above, and the Lie derivative term Lxg in the second component. The formula is well-known 
and actually straightforward to derive, although we note that in |FM2j , for example, the neg- 
ative of the divergence operator is used, so some signs differ; if we consider variations of vr 
as a (l,l)-tensor, then the formula simplifies, as we will note in the next section. ■ 

We now indicate the source of the ten-dimensional obstruction space we mentioned in 
the previous section. 

Lemma 2.4 Let Q be an open domain in with vacuum initial data {6, 0). Then the kernel 
K of the operator is the direct sum of the span Kq of the functions 1, x* (i = 1,2,3) 

and the span Ki of the vector fields Xi = (i = 1,2,3) and = x^-^ — (i < j, 

k ^ i,j)- 

Proof: By Lemma E31 I^$^5,o)(/' ^) = (^1/, Moreover, by Eq.®, LJ/ = 
implies Hesssf = 0, i.e. f G Kq. That LxS = means X is a Killing field of the flat metric 
on M^, so X G Ki. m 

We note that the kernel of D^*^^^^ has a natural interpretation. We first recall that it is 
straightforward to show [Hj that a nontrivial element / in the kernel of L* yields a warped 
product metric —f'^dt^ + g that is Einstein, and if R{g) = then this product metric is a 
solution to the Einstein vacuum equations which is manifestly static. Moncrief [M] showed 
that a nontrivial element (/, X) in the the kernel of D^*^g at a solution of the constraints 
corresponds to a Killing field (symmetry) in the resulting vacuum spacetime. 

3 Constructing Solutions with Good Asymptotics. 

In this section we show that any solution of the vacuum constraint equations with gij = 
6ij + 0(l/r) and vr*-' = 0(l/r^) may be perturbed by an arbitrarily small amount on any 
compact set to a new solution satisfying the condition (AC). We also show that the ADM 
energy-momentum vector is stable under this approximation in that it is perturbed by an 
arbitrarily small amount. This construction shows that the main theorem of the paper ap- 
plies to a very large class of solutions of the vacuum constraint equations which are dense in 
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a suitable sense. 

The construction we are going to make may be viewed as a generalization of the defor- 
mation result of JSYj from the time-symmetric case to the case of arbitrary data. We begin 
with any asymptotically flat data [g, tt) satisfying the vacuum constraint equations and such 
that Qij = 6ij + 0(l/r) and tt = 0(l/r^). We are going to show how to approximate this 
data by new data {g, tt) satisfying the vacuum constraint equations and such that outside a 
compact set we have 

g.. = u'^S-.^ Tc-j = U^{CX)ij (7) 

where u tends to 1 at infinity and X is a vector field tending to at infinity, and where C 
is the operator related to the Lie derivative Lx5 

CX = Lx6 - divs{X)S. 

If such asymptotics can be achieved, then the constraint equations near infinity become the 
equations (computed with respect to 6) 

8Au = u{- |/:X|2 + l/2{Tr{CX)y), AX' + Au-\j{CX)ji - 2u'\,Tr{CX) = 

where the second equation is written with respect to a euclidean basis. Standard asymptotics 
(see 0) then imply that 

u{x) = 1 + a/r + 0{l/r^), X' = b'/r + 0{l/r^) 

for constants a,b\ These asymptotics clearly imply (AC). 

Theorem 1 Let {g, vr) be a vacuum initial data set as above. Given any e > and any 
radius R >> 1, there exists a vacuum initial data set {g,Tf) satisfying ^ which is within 
e of {g, tt) on the ball of radius R. Moreover, the mass and linear momentum of {g, vf) are 
within t of those of {g, vr). 

In order to show that (|7j) can be achieved, we begin by modifying {g, vr) to {g, vr) in the 
annular region from R to 2R so that g = S and vr = outside B2R. We then attempt 
to reimpose the constraint equations by constructing a solution of the form g = u^g and 
vf = v?{Ti + CX) where the operator C is computed with respect the metric g. We hope 
to find a solution u which is near 1 and X which is near 0. The constraint equations then 
become the system 

jj^ = u-^[-8Au + u{R-\n + CX\l + l/2{Trg{TT + CX))^)] = 
{div-g{n))i = u-^[{diVg{7f + CX))i + Au'\{7c + CX){ - 2u~^UiTr^{Tx + CX)] = 0. 

We consider the map T{u,X) = {fl,divg{7r)) and observe that the linearization at (1,0) is 
given by 

DT{ri, Y) = {-8Ar] - Afirj - U'^Yrj + {divY)Tr7c, div{CY)i + Ar]j7ri - 2riiTr7t - 2ridiVg7c) 
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where each term is computed with respect to g. In order to solve the equation T{u, X) = 0, 
we wish to use the inverse function theorem. Thus we would need to check that DT = DTi^i^^^-^ 
is an isomorphism between appropriate spaces with bounded inverse (independent of R for 
large enough R). To be precise, we consider DT as an operator from W'^f to W^f_2- The 
weighted norm convention we are using is that the W'^f norm is given by 

\\fh,p,-s= (/ (i^"/i/^'"')v'rf/^.)^/^ 

0<H<k -^^^ 

where p is a function which equals near infinity and a is a multi- index. It is well known 
that DT is a Fredholm operator of index for p > 1 and 6 G (0,1) (see jB]). Thus in 
order to check that DT is an isomorphism, it would suffice to check that the cokernel is 
trivial; i.e. that DT is onto. This seems to be difficult to check in general, so we use a less 
direct approach. We consider enlarging the domain and we look for solutions of the modified 
problem 

^{u'^g + h, u'^{tt + CX) + A;) = 

where h and k are suitably small symmetric (0, 2)-tensors with compact support. If we can 
solve this problem, the desired asymptotics follow as before. 

We will need to use the fact that the operator D$ = is surjective as an operator 

from ly^f X iyl';_^ to W^°f_2 for 5 G (0,1). This resuh was proven by Choquet-Bruhat, 
Fischer, and Marsden in |("FMj in the maximal case, and by Beig and O Murchadha |B()j 
in the general case. We give a direct proof in the next proposition which appears to work 
under weaker asymptotic assumptions than those of |,BOj . although the same basic idea is 
employed. 

Proposition 3.1 The operator D^(^g .^^ is surjective from the domain W'^f x W^f_^ onto 
W^f_2forp> 1 and 5 G (0,1). 

Proof: We will check the surjectivity for the initial data [gjTc), and the result for {g,Tf) will 
follow by a perturbation argument. We first note the range of D^ contains that of DT (a 
closed finite codimension subspace), and hence is of finite codimension in W^f_2. Therefore 
D$ has closed range. If there is an element (^, Z) in the dual space W^^f^^ which annihilates 
the range of D^, then (^, Z) lies in the kernel of the adjoint and hence satisfies the equations: 

C,^-iA0g^,-^RiJ-lmKP'^Zg),p)g,^ + l/2{K,^ZP)■,p = 
-1/2(Z,;,- + Zj.,i) + {divZ)g,, - 2^Kij + 2i{TrK)gij = 0. 

(Note that these are the adjoint equations with K treated as a tensor of type (1, 1) so they 
differ slightly from those given in Lemma 12.31 The equations in this form can be found in 
CFMj for example, and are straightforward to derive. For the application we have here, 
the detailed lower order terms do not play a role.) We may take the trace and rewrite the 
equations eliminating the terms involving and divZ. We then have 



9 



1/2(Z,., + + 2e^ii = 0. 

Note also that taking the trace of the first equation and the divergence of the second, we get 
a system of hnear equations of the form A(^, Z) = B{x)(y^, VZ) + C(a;)(^, Z) = where 
B, C are coefficient matrices. The following estimates are true for any weight r > and for 
any large radius R 

[ \^p'\'p-'d^^g < c [ (|vvei/+OV"'rf/^. 

JM\Bii Jm\Br 

f {\Z\p^)'p-'dp, < C [ {\Lzg\p'^n'p-'dp, (8) 

where {Lzg)ij = Zi.j + Zj-i is the Lie derivative. A proof of these inequalities will be given 
below. Using these inequalities, we can show inductively that (^, Z) must vanish to infinite 
order at infinity; i.e. \{^,Z)\ < Cjvp~^ for any integer > 1. To see this, we know from 
simple initial asymptotics for the equations that ^ and Z are of order p~^. Putting this into 
the equations we have that VV^ is of order and Lzg is of order p~^. Thus we may 
choose r < 2 and conclude that (^, Z) are of order at most p""^ for any r < 2. Note that the 
pointwise bounds follow from the corresponding weighted bounds by applying standard 
elliptic estimates (mean value- type inequalities) on balls of the form B\x\/2{x). Putting this 
information back into the equations, we find that VV^ is of order p~'^~^ and Lzg is of order 
p~r-2^ Thus we may choose r < 3 and improve the decay on {C,,Z). We thus conclude 
inductively that (^, Z) vanishes to infinite order at infinity. 



To show that (^, Z) vanishes identically, we use a standard unique continuation result 
(see Kazdan ^ for a version written in the most convenient form). To see this, we consider 
doing an inversion (Kelvin transform). That is, we let ) denote asymptotically 

euclidean coordinates and we introduce y = \x\~'^x where |x| is the euclidean norm of x. Now 
observe that the metric g = has scalar curvature R = |x|^(— 8A(|a;|~"'^) + Ixl""*^/?) and 

this is of order at infinity. We then observe that the metric g expressed in the y coordinates 
has Lipschitz components near y = 0, and that gij{0) = 6ij. The conformal transformation 
for the Laplace operators is then Am — l/8Ru = \y\^{A{\y\~^u) — l/8R{\y\^^u)). It follows 
that the quantities (^, Z) = {\y\^^^, \y\^^Z) satisfy a linear system of the form 

A(e, Z) = 5(|/)(Ve, VZ) + Ciy)i^, Z) 

with B is bounded near y = Q and C is bounded by a constant times \y\~'^. The unique 
continuation theorem then implies that (^, Z) vanishes identically. 

Proof of (0): We choose a smooth cutoff function C, which takes values between and 1 with 
C(x) = for |x| < -R and C,{x) = 1 for |x| > 2R. The following inequalities for the euclidean 
metric 5 on i?^, and for functions / and vector fields X in the appropriate weighted spaces 
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are standard and can be proved using integration by parts 

/ \fp^\'p-Hx < C [ {\Vf\p'+n'p-'dx 

JM JM 



{\X\p^fp-^dx < C / {\LxS\p'+^fp-''dx. 

M JM 



Thus we may take / = ^ = to obtain 



/ m^\'p-'dx < c [ {\vm\p'^i'p-'dx 

JM JM 

[ {\CZ\p^rp-'dx < C [ {\L^zS\p'^n'p''dx 

JM JM 

where all quantities are taken with respect to the euclidean metric on the end M\Bpi. Using 
the choice of C, and easy manipulations, these inequalities clearly imply 

/ W^P'^'dx < c! ilV^lp^+^Yp-^dx + CR^^-^ [ ^^dx 

Jm\Br JM\Br JB2r\Br 

[ {\Z\p^yp-^dx < c[ {\Lz6\p^+^)^p~^dx + CR^'-^ [ \Z\^dx. 

Jm\Br Jm\Br Jb2r\Br 

We may now prove (jH)) for the euclidean metric by contradiction. For example, to prove the 
second inequality, suppose we have a sequence Za with 



/ {\Z^\p^fp-Hp, = 1 

JM\Br 



but 



/ {\LzAp'^^fp-'dpg 0. 

J M\Br 



We may assume that the Z^ converge weakly and in L^^^ to a limit Z. By the inequality 
above we see that Z 7^ 0, and we must have Lz5 = 0. Therefore Z is a nonzero Killing 
vector field for 6 with Jj^^^^^{\Z\p'^y p~^dv < 00. This is a contradiction which proves the 
second inequality of (jH)) for the euclidean metric. By a similar argument one proves 

/ m'p-'dx <c I {\vi\p^^^fp-'dx, 

J M\Br J M\Br 

and the first inequality of (jH)) follows from applying this with ^ replaced by |V^| and r re- 
placed by r + 1. Using the fact that g approaches the euclidean metric on approach to cxd, 
the inequalities ^ follow. This completes the proof of the surjectivity of D^. ■ 

We are now in a position to complete the proof of Theorem It is well known that 
{h,k) 1—^ $((7 + /i, vr + /c) is a continuously differentiable map from a neighborhood of (0, 0) 
in W"^^ X W^l_g to W^2-5 for p > 3 and 5 G (0,1). We further observe that the trun- 
cated data (^, 71") are arbitrarily close (for R large) to {g, vr) in W'J X W]:l_^. Thus the 
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linearization D$ is also surjective at {g,TT). We may choose a basis Vi,...,Vn for the 
cokernel of DT, and choose (0, 2)-tensors (/ii, fci), . . . , (/iat, /ctv) in W'^f x W^f^g so that 
D^{hs, kg) = Vs for s = 1, . . . , A^. We now perturb the (/i^, kg) to make them of compact 
support. The corresponding images under are close to the K, and hence they still span 
a complementing subspace for the closed subspace Im{DT). We let W2 be the linear span 
of {hi, ki), . . . , {hj^i, k^). We note that the null space U C W"^^ of DT is A^-dimensional 
since the index of DT is 0. We let Wi be a closed complementing subspace to f/, and we 
let W = Wi X W2, so that is a Banach space. We then define the map T on by 
setting T((u,X), {h, k)) = ^{u'^g + h,u'^{7c + CX) + k). We see by construction that DT 
is an isomorphism at ((1, 0), (0, 0)) e W, so we may apply the standard inverse function 
theorem to assert that T is an isomorphism from a fixed (independent of R) neighborhood 
of ((1, 0), (0, 0)) to a fixed neighborhood of $(^, tt). It then follows that this image contains 
(0,0), so we may find u near 1, X near 0, {h,k) near in the appropriate spaces so that 
$((7, 7f) = where g = u^g + h and n = v^ijr + CX) + k. 

We now show that the ADM energy and linear momentum are perturbed by a small 
amount under this perturbation of (5',vr). If E,Pi denote the total energy and linear mo- 
mentum of a solution [g, tt) of the vacuum constraint equations, we need only show that 
E,Pi are continuous with respect to the W^'J x W^l_g norm for 5 > 1/2. For example, to 
see this for the Pi, observe that by the divergence theorem we have for Ri > R 

j) TiijU^dag - j) TTijU^dag = / T^^jd^g 

J J J Bfi-^\Bji 



|x|=_Ri |x|=-R 



where tt^.^ represents the divergence of the vector field tt^O/OxK Using the constraint equa- 
tion, we see that the volume integrand is a sum of Christoffel symbol times components of 
TT. Symbolically we write 

/ \v\\^Wg < ( / mp'^'Tp-'dpgY'^i [ i\n\p^-'rp-'dpgY^\ 

Jm\Br Jm\Br Jm\Br 

Now if 5 > 1/2 we have 

/ {\rr\p^-'rp~'dpg < ( / {\n\p'^^rp-'dpgf'{ I P'^~'dpgf'-'^" 

JM\Bii Jm\Br JM\Br 

where e = (25 — l)pq/{p — q) > 0. Thus we have 



P - <h TTijU^dag 



< cR- 



|x|=_R 

where the constant c depends only on the W'^f x W^f_^ norm of {g — 5, vr). The continuity 
now follows from continuity of the approximating surface integral. This completes the proof 
of Theorem ^ ■ 
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4 Solving the Deformation Problem. 



In this section we prove the following local deformation theorem for the Einstein constraints. 

Theorem 2 Let Q C be a compactly contained C'^'^'^ -domain, and let qq be a (7^+4>c*_ 
metric and ttq a C''^^'"' -symmetric (2,0)-tensor. Suppose that the linearization D$(g(,^7ro) o/ 
the constraint map $ : C^^'^'°'{Vt) x S^^^^"'{Q) C^''^{Vt) x has an injective formal 

L'^-adjoint D^*^^^.^^^ at (5-0, ttq), where we can consider ' ^tli.^) ^L(^)- 

Then for a sufficiently decaying weight p, there is an e > such that for any function 
V e C*--'"(n) and any vector field W G A''=+^'"(n) for which {{y, VT) -$(^0, vto)) G x 

X^^i''^{Vt) with the support of {(v,W) — $((^0)^0)) contained in Vt and with ^{{v,W) — 

'^'(S'O; ttq)) ||^fc,Q xA''°+^'" ^ ^' there is a C^^"^'"" -metric g and a C''^'^'"' -symmetric tensor n on 

M with ^{g, tt) = (t>, W) in Q and {g, it) = {g^, ttq) outside fi. Moreover, {g, tt) depends 
continuously on {v,W). 

If in addition ((f , W) — ^{go, ttq)) G C^(fi), (go, ttq) and dVl are smooth, and we use an 
exponential weight, then we can solve for {g, vr) smooth. 

If the adjoint of the linearization has nontrivial kernel, then the analogous theorem holds 
for solving ^{g, vr) = {v, W) up to a finite- dimensional error. 

The analogous result for the scalar curvature operator is found in [C]; in particular it 
is shown there that at generic metrics go, for functions S so that the difference S — R{go) 
vanishes outside fl and is sufficiently small (in a weighted Holder norm), there is a metric g 
with R{g) = S, so that g — go is small and supported in f2. The regularity statements are 
analogous to those above. 

The proof of Theorem |21 will be carried out over the next few sections. For simplicity we 
carry out the proof explicitly in the case we require to prove the main theorem, namely we 
have (f , W) = and we are solving up to a finite-dimensional error, as explained below; we 
will then remark on the straightforward modifications for the general case. 

4.1 The Basic Estimate. 

We prove an elliptic estimate in the p-weighted Sobolev spaces in a C^-domain Q. Since 
p decays at dQ, the boundary decay is imposed by p~^, not p, in the sense that tensors 
in p~^-weighted spaces will have to decay suitably at the boundary. We remark that we 
get a global estimate, even though the functions in the p-weighted spaces may not decay 
at the boundary. The key to achieving this estimate is using the overdetermined-ellipticity 
correctly. We state the estimate in the form needed for the asymptotic gluing we will do, 
and later remark how to get the estimate required for Theorem |21 

We let ( G C^{Q) be a bump function, one on most of Q. Let K^, = (K (recall that K is 
defined in Lemma 1231). We define Sg to be the L^((ipg)-orthogonal complement of K^,, and 
we take = Ss- 
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Theorem 3 There is a constant C and an e > so that for all data {g, vr) within e of the 
Minkowski data (5,0) and for all {f,X) G Sg, 

\\{f,X)\\ < C\\D^*(^g^^-^{f,X)\\L2(_n^afig)- (9) 

where the weight p is or e~^^'^^ near the boundary dQ (N > 5/ 

Proof: We first prove estimates on the vector field X. By Lemma l2.ll it suffices to prove 
these estimates for X a smooth vector field on VL. Furthermore we define CgX = Lxg, to 
emphasize we are thinking of the operator acting on X. 

Lemma 4.1 For any weight function p which is d^ or t~^l'^^ (with N >5) near the bound- 
ary dQ, there is a constant C so that for all X orthogonal (or merely in some fixed subspace 
transverse) to the Killing fields of a metric g, the following estimate holds: 

\\X\\hI{Q.) < C\\CgX\\L2(^^y (10) 

Proof: We first note that since the Killing fields are Jacobi fields along geodesies, they are 
locally determined by their first order jet at a point, and so the space of Killing fields on Q 
is finite-dimensional and the fields in the kernel will be smooth (as much as the smoothness 
of the metric and domain will allow) up to the boundary of fl. The Lie derivative is an 
overdetermined-elliptic operator; that is, it has injective symbol. By standard elliptic theory, 
then, one gets an interior estimate of the form 

where Qq CC Q and C depends on d{QQ, dQ). To get a global estimate we foliate Q near dQ 
by level sets Hr = {x ^ Q : d{x, dQ) = r} of the distance function to dQ; there is an tq > 
depending on Q for which these level sets are regular hypersurfaces. By standard elliptic 
theory on these closed hypersurfaces S we get 

II-^eIIhI(S) < C (||£gX||2,2(s) + ||X||i2(2)) 

since the difference between {CgX)\Tj] and Cg^Xj] is zero*^-order in X; gj] is the induced 
metric, and Xy. is the projection of X onto TS. The constant C in this estimate can be 
taken uniform in the distance to dQ sufficiently small. Recall that since |Vc?(-, dQ) \ = 1, the 

r2 

co-area formula gives for positive functions or L^-functions G: j G dpg = j j G dag dr. 

Applying that here yields 

\\iX'^,\/YX'^)\\L2(^A) < G (||£gX||i2(^) + ||X||i^2(^)) 

where is the vector ffeld on A = {x ^ Q : ri < d{x, dQ) < given by the projections 
Xy,. Since the constant is uniform in ri, we see that combining with the interior estimate 
we get 

||(X'^, VsX^)IU2(f70 < G {\\CgX\\L2(n^) + \\X\\L2(n^)) 

where G is uniform in e small, and = {x & Q : d{x,dQ) > e}. We can use the co-area 
formula to integrate the square of this estimate against p'{e) (which is positive for e small). 
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and integrate by parts to get the desired weighted estimate, as in [C]. 



It remains to estimate the terms involving components of X and VX normal to S. We 
compute in a local orthonormal frame 61,62, 63 adapted to S; we take 63 = Vd, well-defined 
in a neighborhood of dQ. We have estimates for Xi-j for i,j = 1,2, and hence by taking 
the trace of CgX, we can also estimate X3.3 as above. It suffices then to estimate Xj.3 for 
i = 1,2. Let Ei denote a set of the form {x G : < d{x, dQ) < rj}, and let C be a cut-off 
function of the distance to the boundary, which is identically one for d < r2 < ri < 1, and 
vanishes outside Ei. We first note 

2 2 

2 / E^'spc^/x, < 2 f cY^xydpg 

I i=i I i=i 

< J ICgX]"^ pdflg - 2 J C'^Xi.^X^.ipdflg. 

ri 771 i=l 



We now want to integrate by parts to make the latter integrand Xi-^iX^; we note that 

by the arithmetic-geometric mean inequality (AM-GM) ab < ^ ^e^a^ -|- ^ j , terms in the 

integrand of the form Xi.^Xj can be replaced (up to a constant factor) by with the 

other term being absorbed on the left side of the above inequality. This observation allows 
us to integrate by parts and switch between covariant and partial derivatives at will. Since 
for z = 1, 2, = = p^j, and since the boundary integrals vanish, we get 



E2 



2 2 

J J2Xl,3Pdf^9 < J \CgX\^pdpg + 2 J CJ2^-'3^X3pdpg + C\\X\\^ 



in)- 



E2 n n 



We now use the Ricci formula to commute covariant derivatives: Xi-j^ — Xi-^j = XiR^j-. So 
we now have 

J^XlspdPg < J \CgX\^pdpg + 2 J CE^v3^3PC?P5 + C||X||i.(^). 

E2 n n 

If we integrate by parts again, we get 

j ^Xl^pdpg < J \CgX\'^pdpg-2 j C^Xi-iXs.^pdpg 

r ^ 

-2 j{Cp),zY,Xi,iX3pdpg + c\\x 



T\\2 I r^w YW"^ 

Hi(n) + ^\\^\\Li{Qy 

1=1 
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We can estimate the second and third summands on the right side above by the AM-GM 
inequahty, bounding them by 

c J |£,x|V/^, + c||x||i.(^) + J \xfd-'pdii,. 

We now claim there is a constant C (depending on N, as well as g and Jl) so that 

X\'^d-'^pdiig <C J \CgX\'^pdng. 

If this is not the case, we can find a sequence X}. of vector fields so that 



/ 



Xk^^d "^pdpg = 1 



but for which 

J \CgXk\^pdpg 0. 

CI 

Now by the previous estimates, then, 

J \VgXk\'^pdpg ij \CgXk\^pdpg + ^ \Xk\'^pdpg I + 1. 

n Vn n / 

Thus we see that X^p^l"^ is bounded in and so we can assume (by taking a subse- 

quence) that there is a vector field X so that Xy^p^l'^ converges weakly in and strongly 
in to Xp^/"^. So X is in Hl^^{Q), and X^ converges to X locally in L^, so that CgX = 
weakly. But the X^ are transverse to the Killing fields, and so X = 0, and thus X^p^^"^ —>■ 
in L^. Now we note that there is an ri and a constant C so that on E^, 

Agp > CN^d-^p. 

So we get by integration by parts 

J C\Xk\^d-^pdpg < 2 J C\Xk\\VgXl,\Nd-'pdpg + ^ C\Xk\^d-'p. (ll) 



Since (' is compactly supported in ^2, the second integral on the right side goes to zero in k. 
By Cauchy-Schwarz we also have 



2 J C\Xk\\VgXk\Nd-'pdpg < 2iV y C\VgXk\'pdpg ^ QX d-^ pd ^g 
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Now using the previous estimate (fTT|). we get that 

j C\Xk?d-^pdiig < 2Ar(l + o(l)) + o(l) 

where o(l) denotes a function which goes to zero as A; ^ oo. So for at least 5, we have 
for large k 

But we also have that 

j{i-0\Xk?d-'p 

n 

goes to zero in k. So far large enough k, 

j \Xk\^d~^pdfXg < 1 
n 

which is a contradiction. ■ 

We continue the proof of the basic estimate. As we noted above f Lemma 12. 4|) . the kernel 
of D^*gQ^ is Kq © Ki. {f,X) G 5*0 precisely if / is orthogonal to (Kq and X is orthogonal 
to (Ki . By the preceding lemma, then, to prove the Basic Estimate at the Minkowski data, 
we now just need to estimate / by L}f. We do have an estimate of the form 

||/lk^(c)<C7(||L;/|U.(o) + ||/|U.(f,)). 

Indeed using Eq.® and taking a trace, we see that 

HessJ = -i (Tr,(L*/) + fR{g)) g - fmc{g). 

Since / is orthogonal to C^o, we can apply standard theory (basically the Rellich theorem) 
to get the estimate 

11/11^.(00 < C||L*/|U.(oo (12) 
where C is uniform in e small, and so we get the desired weighted estimate as above. 
We now have the estimate for {f,X) G Sq 

For data ((?, vr) with e of the Minkowski data, the estimate follows by perturbation. ■ 
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4.2 Variational Method. 



In this section we show how solutions to Ils^D^(^g^T,){h, u) = (0, V) for (0, V) G L^_i(r2) HSq, 
(and {g, vr) close to Minkowski data), can be obtained from standard variational arguments. 

We let denote the orthogonal projection onto 5*0 with respect to the metric S, and 
we define the operator TZig^n) to be the linearization of the map o $, so that 'R'(g,n) = 
Usq o D^(^g .„y We want to show this map surjects from a suitable space onto 5*0 H L^_i(r2). 
In fact it is easy to see this will be true if we can show the map V^g^^) '■= ° D^^g^^^^ is 
surjective to Sg fl L^_i(f2), for {g,^) near the Minkowski data. 

We define the formal adjoint V^g^,) with respect to the metric g, and for p{f,X) G 
L^_i(fi) n Sg we define the functional Q by 



giu,z) = J Q|P(;,)(n,z) r-(/,x) iu,z)^pdfig 



(14) 



We consider the infimum fi = f^(f,x) < over all {u, Z) E Vg = H'^'^{Q) fl Sg. For any 
(m, Z) G Vg, and any {ip, W) with compact support in Q, we have 

{Vl^^iu,Z),{^,W))LH,,^^ = {iu,Z),V^g,^0,W))LH,,,) 

= {{u,Z),D^^g,^){^,W))L2id,.,) 

= (D$^^,,)(w,Z),(^,l^))^.(,^^) (15) 

where the second inequality follows from {u, Z) G Vg. So on Vg we have Z)$^^^^ = V^^^^ 
and thus the Basic Estimate will yield a minimizer to our variational problem, as we now 
recall. We first note that p is finite. 

Lemma 4.2 For any {f,X) G L'^pi^), p > —oo. 

Proof: We simply note that the Basic Estimate, along with Cauchy-Schwarz, yields the 
estimate 

Giu,Z) > ^\\(u,Z)\\l,,,^^^ - ||(/,X)|U.(f,)||(n,Z)||^.a(^). ■ 

Corollary 4.3 For any {f,X) G LKfl), /x = lim Q{ui, Zi) for some sequence {uj, Zj} with 
{\\{.u^,Zi)\\^■2^^^^} hounded. ■ 

Standard functional analysis now gives the existence of a minimizer (mq, Zq) in the Hilbert 
space In fact by the convexity of the functional Q the minimizer is unique. In any 

case, since the orthogonal complement of Sg is composed of C^(n)-functions, we see that 
the Euler-Lagrange equations will hold on all of C^(f2), yielding the weak formulation of 

P(,,.)pP(;,)(wo,^o) =n5,/^$(,,.)p/^$f,,.)(no,Zo) =p(/,X). 
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A simple argument comparing the L^-projections shows that pD^*^^ 7r)('^05 Zq) is also weak 
solution of the equation 

where (0,^) G So fl and where we take P(g^^)pZ}$^^ ^^(mq, ^o) = pi.fi X) to be the 

L^((i/ip)-projection of (0, V") to Sg. 

4.3 Pointwise estimates and the nonlinear problem. 

Assuming we start with tensors {qoiTIq) on M^, sufficiently smooth and sufficiently close to 
Minkowski data on VL CC M^, and which solve the constraint equations outside VL (and in a 
neighborhood of the boundary), we produce tensors vf) which agree with the original data 
outside VL and which satisfy liso^ij), vr) = 0. The basic procedure is straightforward: we use 
the variational procedure to solve the problem at the linear level, and we iterate the process 
of correction. In this section we write down the estimates needed to show convergence of the 
iteration. 

We start by solving the equation IisoD^(^g^^^^)pD^*^g^ .^^-^{uQ, Zq) = -Ilsa^{go,Tro) as 

above. Assuming (go, no) E M'^+^'^'^xS^^^fiJl), we get regularity on the solution [uq, Zq) 
and interior elliptic estimates. We could explicitly derive these, but we will invoke the 
weighting method for elliptic systems of mixed orders due to Doughs and Nirenberg |DJNj . 
Note that for notational simplicity we will omit the subscript [qq, ttq) on the operators below. 

Now the operator p~^D^pD^* is uniformly elliptic, and gives a 4 x 4 system, which 
we will symbolically write as LjU = Ljk{D)U^, and the weights Sj and tk are defined so 
that the order of Lj^ is Sj + tk'-, the elliptic estimate comes from bounding the Holder semi- 

4 4 

norm ^ \U\. „ up to a constant factor by ^ [LiU]_^, ^. In our case U stands for (m, Z), 

i=l " 1=1 

and so we want the weights to satisfy the following conditions: si + ti = 4; for j. A; > 1, 
si + tk = 3 = Sj + ti and Sj + tk = 2. With this in mind we let ti = 4, and the other tk = 3, 
we let Si = and the other sj = —1. Then for Q' GG fl 

\\iho,UJo)\\k+2,a,n' = \\pD^*{Uo, Zo)\\k+2,a,n' < C (||%|U+4,a,n' + ||-^o|U+3,a,n') 

< c(\\iuo,Zo)\\L2(^n) + \\DnipD^''{uo,Zo))\\k,a,n 

+ \\Ddiv{pD<^*{uo,Zo))\\k+i,^^n)- 

If the weight p is near the boundary, then the interior estimates local to the boundary 
(B' = B{x, ^)(ZB = B{x, ^)) are 

\\{ho,OOo)\\k+2,a,B' < Cd^-^'^'^'") {\\{Uo,Zo)\\lHB) + \\p~' Dn{pD<!>* [uq, Zomk,a,B 

+ \\p-^Ddiv{pD^*{uo,Zo))\\k+l,a,B) . 

Here ipi^, «) is linear in k and a, and also note that the lower-order coefficients in the operator 
p~^D^pD^* have powers of 1/d that are accommodated by scaling in the estimate, or can 
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also be subsumed into 'ip{k^a). Using the estimate 

N 

||(mo, ^o)||l2(b) < Cd'~\\uo, Zo)\\L],{n) 

we then have 

+ \\p-^ Ddw{pD(^*{uo, Zo))\\k+i,a,By 

In fact we need to do the above estimates in terms of the projected operator. Because 
the difference between the operators p~^D^pD^* and p~^TZpD<^* on Vg (see fll5|l ) is finite- 
dimensional, we get the following interior estimate 

II (^0, l^o) |U+2,o,f7' = ||P-D$*(Mo, ^o) ||A;+2,a,n' < C* ( || ^0 |U+4,o,f7' + || ^0 ||fc+3,a,C') 

< C (^\{uq, ^o)IU2(n) + \\npD^*{uQ, ^o)|lcfe^- (f^)xA'^^+/.-(c)^ • (16) 

Here Vt' is taken to be a "large" fixed domain containing the support of the fields in K^,. We 
have used that there is a constant C so that for all (0, V) G K^, 

11(0, '^^)llcfe.-(n)xA"=+i.-(n) < C\\{(I),V)\\l2(u>) 

in particular for (0, V) of the form nx,-D$P-D$*(M, Z). We then use that by compactness, 
for all e > there is a C(e) so that 

Z)\\ci{n')y.X^(^') < e||(^, ^)||c"=+4>«(0')X'^'=+^'"(!^') + '^(^)ll(^' ^)\\L^(n')- 

Since the fields in C^K are supported away from the boundary, TZ = near the boundary. 
Thus the estimates we have local to the boundary are as in (fT^. except for the appearance 
of the factor d'^~'^'^^'°'^ on the right-hand side, which yields the decay of the deformation 
tensors at dVt. 

Next we note that we can bound {Jiq^ujo) in L^_i by the basic injectivity estimate for 
1^$*, and we can bound the lower order term in the estimate by the variational inequality 
Q{uq, Zq) < 0, obtaining 

||(/io,o;o)||l2 <C||K,^o)||L?<C||<l>((7o,vro)|U2 

This inequality can be inserted into the decay estimates near the boundary, and it can also 
be used now to get a global estimate on {hQ.ujo): 

\\{hQ,UjQ)\\L2 + ||(/lo,C^o)|U+2,a < C||$(^o,Vro)||^fc.. xA^*+i-"- 

p p-1 p-1 

We now iterate the process of linear correction. We linearize only about the initial metric 
and momentum tensor, since the coefficients in the fourth-order system depend on derivatives 
of ((yfo,7ro); for example in computing LgpL*, one differentiates the Ricci-term twice, which 
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involves four derivatives of the metric. The above estimates show a gain of two derivatives 
in the metric deformation tensor, not four. Hence we cannot apply Newton's method to 
produce a solution. However, we can still produce solutions by linear correction, as we now 
describe. We use Taylor's formula in the form 

$((7o + /i, TTo + a;) = $((?o, vro) + D^^g^^^^){h, uj) + 0{\\ {h, u;)||^+2,J (17) 

where 

110(11 (/i,w)||fc+2.a)||cfc.-xA^fc+i.- . 

sup rrrr — = C{go,iTo) (18) 

and the constant can be taken uniformly for data near {go, ttq) and on subdomains of Q. We 
have the following Taylor's theorem for the projected operator on ^2: 

UsMQo + h,no + u;) = Us.^go, ttq) + 7^(3„,^„)(/^, a;) + 0{\\{h,ij)\\l^^J. (19) 
By solving the equation 

nso^-^* (9o,-o)P^'^ (90,-0) ("o,^o) = -H5o$((7o,vro) 
as above, then, we have that 

H5o$(^7i,7ri) = 0(||(/io,^o)|lU«) 

where {ho,uo) = pD^*^g^^^^~^{uo, Zo) and (5-1,711) = {go, no) + (^0,^0); this holds locally near 
the boundary as well, outside the support of (. The quadratic decay at the first step does not 
propagate (since we linearize about the initial data only) but the estimates above show that 
a Picard iteration converges; indeed it is straightforward to check that the error terms we get 
by using the fixed linear operator to correct the nonlinear term at each stage still allow the 
iteration to converge geometrically at a rate better than linear, but worse than quadratic, 
adapting the proof of Prop. 3.9 from (Cj. (We remark that there are a few harmless typos 
in [C] ; in particular, C^l" should be defined as we have done here, and so the norm one uses 
on {h,uj) (simply h in 0) in the Picard iteration is the norm on the left side of Eq. (j2(Jj) 
below.) The local estimates near the boundary (outside the support of () allow us to show 
the limiting tensors can be made to decay as much as we like as we approach dfl by choosing 
large enough; we can even get all derivatives decaying if we choose smooth data ^{go, ttq) 

sufficiently small, supported away from dQ, and if we use an exponential weight p. Moreover 

00 

we have a bound on the limiting tensors {h, uj) = J2 i^k, ^k)'- 



k=0 



\\{h,u;)\y^^ + \\{h (20) 

This bound comes straight from the estimates done above (and the geometric convergence). 

We remark that for the smooth case we solve the nonlinear equation in some finite regular- 
ity class to start, but using an exponentially decaying weight. We use standard bootstrapping 
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on the quasilinear elliptic condition $(((7o,7ro) + pD^*{u, Z)) G -fC* C C°° to see that the 
solution pD^*{u, Z) is smooth on the interior, and then use the decay of p to prove that 
all the derivatives decay near the boundary; we note that we have assumed that ^{go, ttq) is 
supported away from the boundary, and we are considering the case {v, W) = of Theorem 
El The general case is similar, as we discuss in the next remark. 

Remark. We now note what modifications are needed to carry out the above argument to 
prove Theorem |21 as stated above. Essentially we just need to derive the analogous Basic 
Estimate. The variational method and nonlinear iteration will proceed as above. In fact, 
when there is no cokernel present, the method is slightly easier, since we do not need to worry 
about staying transverse to a subspace. The Basic Estimate follows by first observing the 
estimate ||X||j:/i(q) < C(||£c,X||i2(Q) + ||X||^2(f^)) follows by combining the estimate above 
with the fact mentioned earlier that the kernel of Cg is finite-dimensional and smooth up to 
the boundary. Note that by differentiating the momentum constraint (and using the Ricci 
formula) we see that the condition that (/, X) is in the kernel of D^* reduces to a second- 
order ODE system along geodesies, and so the kernel is also finite-dimensional and smooth 
(as smooth as allowed by the other data). Elementary manipulation using Lemma f2.3l then 
yields the estimate 

||(/,X)||^2a(^) < c{\\D^l^){f,X)U.^n) + ||(/,X)|U2(^)). 

We can remove the lower-order term in the estimate by essentially the usual Rellich argu- 
ment; to do this, one first shows the left side of the above estimate can be replaced by 
II (/, X)p2 ||j^2(Q). In fact we already have the required estimate for the vector field X above, 
and the analogous estimate for / in terms of L*f is a straightforward modification; if we 
had incorporated different weights on the derivatives in the definition of the weighted spaces, 
this would have already been included in the Basic Estimate. This completes the proof of 
Theorem |21 ■ 

5 Handling the cokernel. 

We apply the preceding analysis to solve the constraint equations with given model near 
infinity up to a finite-dimensional obstruction. Given any AF solution of the constraints on 
M.^, we fix AF coordinates at infinity, and at a sufficiently large radius R we smoothly patch 
(using a smooth cut-off function) our original data to the model data, the slow-patching oc- 
curring in the annulus Afi from R to 2R. To apply the preceding section, we scale our glued 
data to the unit annulus Ai by ((/^(x), 7r/j(x)) := {g{Rx), R7i{Rx)). Under this scaling the 
relevant geometric quantities (derivatives, curvatures, etc.) are 0{R~^), the (AC) conditions 
hold to order 0{R~^), and the vacuum constraint equations are scale- invariant. For large 
enough R the data will be sufficiently close to the flat data on the fixed annulus Ai, so the 
data can be perturbed to make the constraint functions lie in (K on Ai. 

It is at this point we need to be clear about what are suitable families of solutions to 
glue on near infinity. We define a family of solutions on the exterior of a fixed ball and 
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smoothly parametrized on an open set O C M^^ to be admissible if, with reference to a 
fixed coordinate chart near infinity, the family satisfies (AC) locally uniformly, and the map 
B : (9 — > M.^^ which associates to each member of the family its energy-momenta [E, P, J.C) 
is a homeomorphism onto an open subset. 

We note that slices in Kerr form an admissible family, parametrized by the total mass 
m, the angular momentum parameter a, and an element in the Poincare group to control 
Euclidean motions of the AF coordinate system as well as boosts. It is well-known that 
under the action of the Poincare group the energy-momentum transforms as a four-vector, 
and the center and angular momentum enjoy a similar transformation rule (they comprise 
a skew-symmetric M^,^ with Mq^ a constant times C^) cf. |ADMj . |BUj . |RTj . which allows 
one to see that varying the slices near the given one gives a local homeomorphism from a 
ten-dimensional family of slices to the ten-dimensional space of energy-momenta. Note that 
time translations and rotations about the axis of symmetry of the Kerr are isometrics, so 
we really mod out by the closed two-parameter subgroup they generate to parametrize the 
slices, leaving eight effective parameters from the group. Please see )CDj for an explicit 
formula and other examples. 

Let C M correspond the the exterior {x G : |x| > R} in an AF chart. We now 
state the Main Theorem: 

Theorem 4 Let ((7,7r) solve the vacuum constraint equations on M^, and satisfy (AC), in 
an AF coordinate chart in a given end. Let O parametrize an admissible family of solutions, 
and Aq G O has matching that of {g,n). There is a radius R and a solution {'g,W) of 

the constraints so that (^, vf) = ((7,7r) on M \ Er, and {'g,7f) agrees with a suitably chosen 
member of the admissible family on E2r. 

Coupled with Theorem ^ we have the following approximation result. 

Theorem 5 Let [g, tt) be any AF solution of the vacuum constraints. Given any e > 0, 
there is a solution (^, vf) within e of ((?, tt) and whose ADM energy-momentum {E,P) is 
within e of that of [g, tt), so that near infinity, {g,Tf) agrees with a member of an admissible 
asymptotic model family. 

Proof of Theorem ^ For sufficiently large R, we have a continuous map X from an open 
set O C M}^ to as follows: take data corresponding to A G O, glue it to the given 
data in Ar, and then perturb using the above techniques to data {g,7i) with $(^,7f) lying 
in a fixed ten-dimensional vector space; we identify $(^,7f) with X(A). By design, (^, vr) is 
identical to the member of the asymptotic model family in the exterior end E2R. In the next 
section we show that this map X has a zero near the parameter Aq G (9 corresponding to 
the given initial AF data [g, it) with which we started, and thus this data solves the vacuum 
constraints. Notice that we cannot arbitrarily glue on anything we like, but the procedure 
finds a suitable exterior that will lead to a solution of the constraints. 
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5.1 Computing the parameter map. 

In this section we verify that the compactly supported deformation which puts the constraint 
data into the approximate kernel preserves the asymptotic symmetric conditions (AC). This 
allows us to conclude that the projection of the data onto the cokernel is given to leading 
order by the change in parameters across the annulus from the initial to the model we 
glued on; i.e. the deformation does not induce any extra error terms. This simplifies the 
computations in [Cj, where less precise estimates were used along with explicit calculations 
in the conformally flat case. 

5.1.1 Preservation of the AC Condition 

Starting with g and vr satisfying the asymptotic condition (AC) (in an AF chart, includ- 
ing several derivatives), we glue on a model solution on Ar, using a symmetric cut-off 
function which is zero near one of the boundary spheres, and one near the other; this pro- 
duces an approximate solution (g,TT). We note that the model solutions which we con- 
sider satisfy (AC), uniformly near a given one. Next, we scale to the unit annulus Ai: 
{9R)ij{x) '■= 9ij{Rx), and {nii)ij{x) := Rnij{Rx). We see that {gn)ij{x) - {gR)ij{-x) = 
0{R~^), {iTR)ij{x) + (iTRjiji—x) = 0{R''^), and similarly the derivatives will satisfy the 
appropriate even/odd condition to order 0{R~'^), and hence the Christoffel symbols and 
curvatures do as well. We use the notation u{x) = u^{x) = u{—x), and we also note that 
"O" will include several derivatives of the quantity in question (as required in the next sec- 
tion), and on the annulus Ai the derivatives will decay at the same rate by the scaling. 

Let {h,uj) = pD^*^^^^^-^{u, Z) solve the nonlinear projected problem; in particular, 

oo oo 

{h,u) is obtained by iteration, so that {h,uj) = J2i^k,^^k), where {u, Z) = J2{'^k,Zk) 

k=0 k=0 

and {hk,uJk) = pD^*^g^^^^^{uk, Zk), and 

3 3 3 

"^{gR + h,7rR + Uj) = Ci^^'C, ^^^^^ + ^kYkC) 

i=0 j=l k=l 

on Ai, with x^ := 1. Since (^{gn^TTn) = 0{R-^), we have by ^ that {h,uj) = 0{R-^). 

We now note that the anti-symmetric part of ^{gR + h,7[R + uj) is small. Indeed, for 
i = 1, 2, 3 we have d = 0{R-'^), h = 0{R-^y. 

[ nign + h,nn + uj)x' dx = [ [Rign] + Lg^ih) + Oi\\h\\la) + OiR-^)]x' dx 

JAi JAi 

= / [R{gR) - RM]x' dx + / Lg^{h)x'dfig^ + 0{R-^) 

Ja+ JAi 

where we have used the fact that \dx — dfig^\ = 0{R^^). The first term above is 0{R^'^) by 
the symmetry of R{gR), as is the second integral, since L*^(x*) = 0{R~^); indeed we note 

{Hessg.x^)^^ = ~T]lidx% = OiR-'). 
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Similarly, using the symmetry of Yk 

/ diVgj^j^h{T{R + uj) -Yi dx = I divs{TTji + Lu) ■ Yi dx + 0{R~'^) = 0{R~^). 
Jai Jai 

It directly follows that D<l>(g^^^^)(/i, — -D$(g^_7r^)(/i, x) = 0{R^'^). With a little 
computation, we also see that D^(^gj,^^j,){h-h, uj+uj), p((-D$(^^ Z))^-D^*^g^^^^~^{u, Z)^ . 

and hence D(^(^g^^.„^)pD(^*^^g^ ^^^{u — u, Z — Z), are also 0{R^'^). Here we use the fact that we 
have uniform bounds on pV^u (/c = 0, 1, 2) and pV^Z (A; = 0, 1) and their derivatives from 
the elliptic estimates (and the decay of p). 

Now we estimate {h-h,uj + u)= p(^D^l^^^Ju,Z) - {D^l^^^Ju,Z))''^+0{R-^) = 

p(jJ^*g^ ^^){u, Z) - D^lg^ TTfl)!""' ^)) + 0{R~^). By the elliptic estimate, using the interior 
estimate near the boundary to establish decay and global bounds, we see that it suffices 
to show that {u,Z) — {u, Z) is 0{R~^) in L^(y4i). To show this, we note that the above 
symmetry discussion, coupled with linearization and the decay estimates obtained from the 
elliptic estimate, shows that D$(g^,7r^)pD$^^^ TTji)i'^ — u,Z — Z) G (Ai) and that we can 
write 

/ ((n, Z) - i^))D^^g^,^^)pD^l^^^^Ju -u,Z- Z)dpg, = \\ {u, Z) - K^|U.0(i?-2). 

J Ai 

Furthermore, (m, Z) and hence (m, Z) are transverse to K^. This allows us to apply integration 
by parts (using Lemma f2.ip and the basic elliptic estimate transverse to the cokernel to get 

/ Z) - (m, Z)) D<l>(^g^^^^)pD<!>lg^^^^^ {u-u,Z - Z)dpg^ 

JAi 

= [ p\D%^^^Ju-u,Z-Z)\'dpg, > C\\{u,Z)-'^)\\ 

JA^ 
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5.1.2 Controlling the resulting constraint data. 

The preservation of the (AC) condition implies that the projection of the constraint functions 
in the direction of the cokernel after solving the nonlinear projected problem is given, up 
to lower-order error terms, by boundary integrals; of course the boundary data and its 
derivatives are unchanged by the perturbation in the annulus. In particular, this means that 
the constraint functions can be controlled by varying the parameters of the solution we glue 
on, and so as long as the family we glue on has enough degrees of freedom, we can make the 
Hamiltonian and momentum constraints zero. Indeed, we let {gRiT^R) = {qr + /i, ttj? + a;), 
we let (^, 7f) be the re-scaled version on Ar, and we note for A; = 1, 2, 3 

/ x''n{gji,Wii)dx = / x''y2[i9R)ij,ij - {9R)iLjj]dx + 0{R-^) 
Jai Ja-, 

= R-^ I x^y^ idwj - 9n,n) dx + 0{R-^) 
Ja„ T7 
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Here we have used the fact that the difference between -R(^^) and ^ [{gR)ij,ij ~ i9R)ii,jj\ 

0{R^^) and is even to 0{R^^), and similarly for the quadratic terms in Wr in the Hamiltonian 
constraint. We now note the other projections, starting with the mass term, the projection 
onto the constant functions: 

/ U{gR,TtR)dx = / y^[{gR)ij,i3-{gR)ii,33\dx + 0{R~'^) 

J A, J A, 

= [ J2i9ij,iJ-9u,n)dx + 0(R-') 

Integrating the momentum constraint against the translation fields Xj. yields 
/ diVgj^(WR) ■ Xkdx = / divs(WR) ■ Xk dx + 0(R~^) 

JAi J Ax 

= / {j:R),rXlv' di + 0{R-') 

JdAi 

= R-^ [ WijXy dc + oiR-"^). 

Finally, we project the momentum constraint onto the rotation fields Yk'. 

/ diVg^(WR) • Yfe = / divs(WR) ■Ykdx + 0(R~^) 

JAi J Ax 

= / (7f«k-nv o?e+o(i?-3) 

JdAi 

= R-^ [ WijY^u^ d^ + 0{R-'). 

JdAn 

Note that we could also have used the measures from the metric in these integrals. 
Also recall that (^, 7f) agree with the given data {g, tt) on the inner boundary, and agree with 
the data of the asymptotic model on the outer boundary of the annulus. To show the map 
X has a zero, we show the associated map : A i-^ 

[r I n{gii,nR)dx,R diVg^{nR)-Xk,R^ diVg^{nR) ■ Yk dx, R^ x^n{gR,TfR)dx) 

^ JAi JAi JAi J Ax ' 

has a zero. 
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5.2 Solving the constraints. 

We now finish the proof of the Main Theorem. Thus far, we have not chosen which member 
of the model family {e.g. which slice in which Kerr) to use near infinity. We understand that 
the model will have ADM energy-momentum, and angular and linear momentum, near that 
of the original data. Suppose Aq G parametrizes a solution whose energy- moment a 0(Ao) 
agrees with that of {g,7f). We then consider A near Aq. 

The computation in the previous section shows then that the map Ir{\) = 0(A) — 0(Ao) + 
0(1), where o(l) ^ as i? ^ 00, uniformly for A near Aq. Fix a small ball B about Aq. Let 
2ji(t, A) = 0(A) — 0(Ao) + to(l) be a homotopy defined on [0, 1] x B, between the homeomor- 
phism 0(A) — 0(Ao) and the map Xr. For sufficiently large R, we have ^ Ir{[0, 1] x dB). 
By degree considerations |Nj, for large R, Xr must hit zero for some A G -B. For such a 
A, the constraints are satisfied, with the model near infinity corresponding to A for which 
= 0, and with 0(A) near 0(Ao). ■ 
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